function answer=integral(x1,x2)
div=1000
delta=(x1-x2)/(div-1)
answer=0
for i=linspace(x1,x2,div)
    a=f(i)
    b=f(i+delta)
    s=(a+b)*delta/2
    answer=answer+s
end
